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Abstract 

Some probabilistic aspects of the number variance statistic are investi- 
gated. Infinite systems of independent Brownian motions and symmetric 
a-stable processes are used to construct new examples of processes which 
exhibit both divergent and saturating number variance behaviour. We de- 
rive a general expression for the number variance for the spatial particle 
configurations arising from these systems and this enables us to deduce 
various limiting distribution results for the fiuctuations of the associated 
counting functions. In particular, knowledge of the number variance al- 
lows us to introduce and characterize a novel family of centered, long 
memory Gaussian processes. We obtain fractional Brownian motion as a 
weak limit of these constructed processes. 
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Introduction 

Let (X, JF, P) be a point process on M, that is, a collection X :— {{xi)°^_^ : 
Xi € R, Vi and #(2:; : x., e [6, 6 + i]) < 00 V6 G K, L G M+}, with the min- 
imal cr-algebra generated by these point configurations and P some probability 
measure on {X, T) . The associated number variance is defined as 

Var{N[b,b + L]) := Variancep[#(a;i : Xi € [b,b + L])]. 

More generally, in order to deal with non-spatially homogeneous cases, it is more 
convenient to work with the averaged number variance which we define as 

V[L] ■.^Eb[Var{N[b,b + L])] 
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taking an appropriate uniform average of the number variance over different 
intervals of the same length. As we increase the length L of the underlying 
interval the number of points contained in that interval will also increase. How- 
ever, in many situations, it is not immediately clear what will happen to the 
variance of this number as the interval length grows. One of the main questions 
considered in this paper will be the behaviour of T^[L] as L ^ oo. We shall see 
that, somewhat counter- intuitively, in some instances we have 

lim V[L] — K < oo, 

L — 'OO 

in which case we will say that the number variance saturates to the level k g . 

Until recently, interest in the number variance statistic had been confined 
to the fields of random matrix theory (see e.g. ^lEl) and quantum theory 
(e.g pi I17p where it is commonly used as an indicator of spectral rigidity. In 
the study of quantum spectra, the manner in which the number variance (of 
eigenvalues) grows with interval length provides an indication of whether the 
underlying classical trajectories are integrable, chaotic or mixed. In the large 
energy limit, the spectral statistics of a quantum system with strongly chaotic 
classical dynamics, should, according to the conjecture of Bohigas, Giannoni 
and Schmidt [7], agree with the corresponding random matrix statistics which 
are independent of the specific system. However, in reality, for many such sys- 
tems, when it comes to the long range "global" statistics, this random matrix 
universality breaks down. In these cases, the number variance typifies this tran- 
sition to non-conformity in that following an initial random matrix consistent 
logarithmic growth with interval length, the number variance then saturates. 

Attempts to improve the understanding of the deviations from random ma- 
trix predictions have led to convincing explanations for number variance satu- 
ration behaviour in terms of periodic orbit theory, see for example 0] ■ In 
^ there is a heuristic derivation of an explicit formula for the empirical number 
variance of the zeros of the Riemann zeta function on the critical line which is 
consistent with numerical evidence. 

In the last few years, number variance has been considered from a slightly 
different viewpoint in relation to determinantal point processes (see (24. for 
background on this topic). Results on the growth of the number variance for 
these processes are given in e.g. EHl- The emphasis in these cases is on 
ascertaining the divergence of the number variance as this, it turns out, is the 
key ingredient needed to prove Gaussian fluctuation results for the counting 
functions of a large class of determinantal processes, including those arising 
naturally in the context of random matrix theory. 

Motivated by the Riemann zeta example, Johansson recently introduced 
an example of a point process with determinantal structure which demonstrates 
the same type of number variance saturation behaviour as conjectured by Berry 
(El) for ftic Riemann zeroes. This process is constructed from a system of n non- 
colliding Brownian particles started from equidistant points Uj = T + a{n — j) 
with T £ R, a G E+ j = 1, . . . ,n. There are a number of approaches to 
describing such a system, see El for details. In any case, it can be shown 
that the configuration of particles in space formed by the process at a fixed 
time is a determinantal process and as such its correlation functions (or joint 
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intensities) take the form 



i?(;;)(xi,a;2,...,x™) = clet(ifl"^(a;„a;j))" . (0.1) 

\ / i,j=i 

Here and for determinantal processes in general, the determinant may be in- 
terpreted as the density of the probabihty of finding m of the points in the 
infinitesimal intervals around Xi,X2, ■ ■ ■ , Xm- 

As the number of particles n — > cx), Johansson shows that the correlation 
kernel K^^'^ converges uniformly to a kernel with leading term (the correction 
is of order e~'^^,(3 > 0) 

1 /smTT{x-y) dcosn{x + y) + {y-x)smn{y + x)\ 

Kt(ax,ay} = -[ — + — , (0.2) 

a V 7r(x — y) 7r(d^ + [y — xj'^j / 

where d = This limiting kernel defines, at each time t, a limiting, non- 

spatially homogeneous determinantal process which may loosely be thought of 
as the configuration of particles in space formed by an infinite system of non- 
colliding Brownian particles started from an equispaced initial configuration. 
When the interval length L is small relative to d, the averaged number variance 
for this process computed from the leading part of the correlation kernel (|0.2(l 
has leading term 

^(log(27rL/a)+7E„,„ + l). (0.3) 

However, if d is held constant, while L is increased, it is deduced that the number 
variance saturates to the level 

^(log(27rd)+7,„,„ + l). (0.4) 

The "small L" expression ()().3|l agrees with the number variance of the determi- 
nantal point process associated with the sine kernel of density a 

i^(..,.,-)- """/"-'"f^ (0.5) 
n[Xi - Xj) 

This determinantal process is the universal scaling limit obtained for the eigen- 
values of random matrices from e.g. the Gaussian Unitary Ensemble and U(n), 
as matrix size tends to infinity (see for example, |ll9p. 

For our purposes it will be convenient to think of the above averaged number 
variance as the number variance of the spatial particle configurations arising 
from an "averaged model" which we choose to interpret as an infinite system of 
non-colliding Brownian motions started from the initial positions 

Uj^a{j-e), j el, ae]R+, e - Uniform [0, 1]. (0.6) 

The two key features of Johansson's model as just described are the equis- 
paced starting positions and the fact that the Brownian particles are conditioned 
to not collide. In this work we consider the number variance statistic in the in- 
dependent process analogues of Johansson's model. Since these independent 
versions do not have determinantal structure, existing number variance results 
no longer apply. 
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The paper is organized as follows. We begin by deriving a general formula for 
the number variance for the spatial particle configurations arising from infinite 
systems of independent Brownian motions and symmetric a-stable processes in 
K started from the initial configuration 1)0.6(1 . This enables us to deduce the 
asymptotic behaviour of the statistic as the length of the interval over which 
it is defined goes to infinity. We give an explicit formula for the saturation 
level for the cases in which saturation does occur. Once this is achieved we 
are able to explain the number variance saturation phenomenon in terms of 
the tail distribution behaviours of the underlying processes. We provide two 
specific illustrative examples as corollaries and briefly discuss how our results 
compare to those already existing in the literature. We conclude the first section 
by demonstrating the close relationship between the number variance and the 
rate of convergence of the distribution of the associated counting function to a 
Poisson law. 

In the second section we use the number variance to prove Gaussian fluc- 
tuation results for the counting functions of our particle conflgurations in two 
different scalings. In the third and flnal section we add some dynamics to the 
fluctuations of the counting functions to construct a collection of processes, each 
of which is shown to converge weakly in C[0, oo) to a centered Gaussian pro- 
cess with covariance structure similar in form to that of a fractional Brownian 
motion. Our earlier results on the behaviour of the number variance allow us 
to better characterize these limiting processes. In particular, the long-range 
dependence property exhibited by the covariance of their increments is directly 
determined by the rate of growth of the associated number variance. In the 
cases corresponding to a G (0, 1), a further rescaling of the limiting Gaussian 
processes allow us to recover fractional Brownian motions of Hurst parameters 
as weak limits. 

1 The independent particle cases 

1.1 A Poisson process initial configuration 

We begin by illustrating the effect of the initial positions on the number vari- 
ance of the spatial particle configurations arising from such infinite systems of 
processes as those considered in this paper. The following theorem is the well 
known result (see for example, Chapter VIII, section 5) that the Poisson 
process is invariant for an infinite system of independent particles with the same 
evolution. 

Theorem 1. Consider an infinite collection of independent, identical in law, 
translation invariant real-valued stochastic processes {{Xj{t),t > 0);j G Z}. 
Suppose that {Xj{0)}°^_^ is a Poisson process of intensity 9 on M. Then 
{Xj(t)}°°^_^^ is a Poisson process of intensity 9 for every t. 

Consequently, if we begin with a "mixed up" Poisson process initial configu- 
ration and allow each particle to move as (say) a Levy process, independently of 
the others then we observe a non-saturating, linearly growing number variance 
{V[L] — 9L) at the start and for all time. 
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1.2 The Brownian and symmetric a-stable cases 

This last theorem served to highhght the importance of the regularity or rigid- 
ity of the starting configuration H0.6|l in determining the number variance be- 
haviour. However, it is reasonable to suppose, that in Johansson's model, the 
strong restrictions placed on the movement of the Brownian particles must also 
contribute significantly to the saturation behaviour that its number variance 
demonstrates. This leads us to ask; what would happen if we started with an 
initial configuration such as (|().6|l but did not place such strong restrictions on 
the movement of the particles? We answer this question for the cases in which 
each particle moves independently as a one-dimensional Brownian motion or as 
a symmetric a-stable process on M. 

Recall (see e.g. that a symmetric a-stable process is a Levy process 

{X"''^{t),t > 0) with characteristic function, for each t>0, given by 

<Ax=.c(t)(0) E[e^«^°-^W] 

= exp(-tc|6'|"), oO, ae(0,2). (1.1) 

Some of the properties enjoyed by this class of processes are; 

• {X"''^(t),t > 0} with associated transition density pt{x,y), x,y € M. is 
temporally and spatially homogeneous. 

• {X°'''^{t)} is symmetric and self-similar in the sense that {X"''^(t)} 
{-X"''=(i)} and {A-i/"X"''=(At)} {X"''=(i)} for constants A G K. 

The arguments that follow also apply to the a — 2 Gaussian cases. Note that 
we have standard Brownian motion when a = 2, c = ^. 

Theorem 2. Fix a symmetric a-stable process (X°'''^(t), t > 0) on R with prop- 
erties as described above. Suppose we start an independent copy of this process 
from each of the starting positions 

Uj := a{j — e), j € Z, where a £ and e Uniform[0, 1]. (1-2) 

The configuration of particles in space formed by this infinite system of inde- 
pendent symmetric a-stable processes at a fixed time t has number variance 

Proof. Let {{X°''''{t),t > 0),j G Z} be the independent copies of the chosen 
symmetric a-stable process indexed by j € Z. Denote the law of each X"''^ 
started from x G M by P^, and write P := P*^. Now the number of symmetric 
a-stable particles in an interval [0,L] C M at time t is given by the sum of 
indicator random variables 

oo 

N^'^^L]^ ^ I[Xf'^{t)+u,e[0,L]], (1.4) 

j = -oo 

where {uj)°°^_^ is given by (|1.2(l . Note that by construction, for this "averaged 
model", for all 6 € R we have 7V"''='"[0, L] 7V"''='"[6, 6 + i]. Thus the number 
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variance is given by 

14"'"'" [L] := Var[A^f'"'°[0,L]] 

OO 

= Y.^[X';'\t)+u,e[0,L]]nX';-'{t)+u,i[0,L]]. (1.5) 

j = -oo 

We can use the self-similarity property and the independence of e and the Xj 
to write the probabilities under consideration as convolutions which then allows 
us to deduce 

^ P[X;'=(t) + a{j -e)e [0, L]]= / pt/,. (x, y) dy dx. 

Hence 

y.-'^'^L] (1.6) 

L/a poo pL/a pL/a 

Pt/a''{x,y)dydx~ / / / Pt/a''{x,y)pt/a''ix,z)dzdydx . 

J-ocJo Jo 



Ti T2 

By Fubini's Theorem and symmetry we have 

f*L/a /"OO 



Ti= / pt/a''{y,x)dxdy ^ -. 

Jo J-oo a 



(1.7) 



For the other term we make use of the Chapman-Kolmogorov identity and the 
spatial homogeneity before performing an integral switch to obtain 

rL/a p{L/a)~y 

T2= I I P2t/a"(0,^)dzdy 

Jo J-y 

pL/a I'O rL/a 

= - / P2t/a°(0,z)dz+ / Zp2t/a"(0,z)dz- / Zp2t/a° (0, z) dz. 

J-L/a J-L/a Jo 

From the symmetry property we know that, for each pt(0, z) is an even func- 
tion in z and g{z) := zpt(0, z) is an odd function in z. These facts allow us to 
conclude, bringing the two terms together. 



'[L] = / P2t/a°(0, z)dz + 2 / Zp2t/a°(0,2:)dz 

Jhla Jo 



(1.8) 



Applying Fourier inversion to the characteristic function 0x° ':(t) (^) given at 
()1.1|) . we deduce that the transition density can be expressed as 



Pt(0, zja, c) = — / cos(z(?)e 



dQ. 



Using this density, the symmetry property and the expression 1)1. 8|l we obtain 



\L\ - - + - 
a IT 







OO p-2ct(e/a)° r L9. 

cos ( — ) — 1 

a ' 



de. 



as required. 



□ 
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Having found a general expression for the number variance we are able to 
consider its behaviour as the interval length L is increased. Note that given 
positive real valued functions g and h we lei g ^ h signify that lim ^ — 1. 

Theorem 3. Consider the number variance V^'''^'°'[L\ for the system of sym- 
metric a-stahle processes considered above. We claim that 



ka,c^ajL^~'', for a e {0,1), 

ka,c,a,t'^og{L), for a ^ I, 

ka,c,a,t + Ksat{a, c, a, t) , for a G (1, 2), 

ka,c,a,t Z^"*- + Ksat{a, c, a, t), for a = 2, 



(1.9) 



as L ^ oo, where ka^c,a,t is a constant depending on a, c, a and t, and 
Ksat{a,c,a,t) = — (2te)i/"rfl - -), 



a 



[1.10) 



where T{x) :— s^ ^ e ds is the usual Gamma function. 



Proof. From (|1.8fl note that we may re-write the expression for the number 
variance as 



L 

a 



X..(5)|>i/<.]+E,[|x-(5)|;|x«(f)|<L/ 



(1.11) 



Now the behaviour of the first term in this expression is well known (see |2()| . 
page 16). For a e (0, 2) we have 



L 
a 



a 



l-a 



as L ^ oo. 



where ka — x " sinxdx) . We use log(L) when a — 1. In the Gaussian 
case (a = 2) we have instead 



L 
a 



X 



2t\ 



> L/a 



V2c7rf 



:e-^/8^* as L 



To deal with the second term of observe that we have 



X"'"(2i/a") ; X"'%2t/a°') < L/a 



X°''''(2t/a") 
X°'^'={2t/a°') 



> X 



dA. 



Thus it is clear that the rate of divergence or saturation of the number variance 
is determined by the upper tail of the distribution of the underlying symmetric 
a-stable process. This may also be seen by differentiating the integral expression 
for Vt"''^'''[L] with respect to L 



dL 



> 
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Consequently, for a G (1,2] (the saturating cases) 



lim V° 

L — *-oo 



X 



Ksat{a,c,a,t) < oo. 



The exact expression for Ksat can be obtained from the moments. By j22j . if X 
is a symmetric a-stable random variable with < a < 2 and scale cr, then for 
— 1 < ^ < a we have 



E 



X 



^s/a 2'5r(i±^)r(l-(5/a) 



r(i/2)r(i-<5/2) 



Applying this theorem with 5 = l,a 
with the integral expression for V^"'"^'' 

Ksatia,c,a,t) = - 



= ^ gives (jl.lUII . To see how this fits in 
'[L], it may be verified that for a e (1,2] 

2 poo g-2ct(e/a)°' 



■ de. 



□ 



The saturation result is now a consequence of Theorem |21 

Corollary 1. 1. If a E (0,1], at each time t > 0, V"''''"'[L] diverges 
L — > oo. 



2. If a e (1,2], at each time t < oo, 
i^sat{ct, c,a,t) as L — > oo. 



'[L] saturates to the level 



Remark 1. Even from the simplest equation H1.5|l it is clear that the largest 
contributions to the number variance come from the activity at the edges of 
the interval under consideration. Thus intuitively, the fatter the tails of the 
distributions concerned the greater the number of particles that may be in the 
vicinity of these edges making these substantial contributions, the slower the 
decline in this number as the interval length increases and consequently the 
slower the decay in the growth of the number variance. 

We now apply Theorems |2] and |2| to two well known examples. 

Corollary 2 (Brownian case). Consider an infinite system of Brownian par- 
ticles started from the initial configuration {uj)°°^_^ as given at H1.2|l . The 
number variance for this process is 



a 



1 - e-^'/4* 



whe 



1 



dy. 



J-oo 

As L ~* oo this number variance saturates exponentially quickly to the level 

t 



Corollary 3 (Symmetric Cauchy case). Consider an infinite system of sym- 
metric Cauchy processes started from the initial configuration {uj)^_^ as given 
at (|1.2I) . The number variance for this process is 

V,'-'-[L] = ^ [l - ^ arctan [^]\ + ^ [log (l + . (1.12) 

a I TT 2t i aiT I V \2t/ / i 

For "large" (relative to 2t) L, we have 



OTT \2,t/ 

and so the number variance diverges at a logarithmic rate as L ^ oo. 

Remark 2. Coincidentally, in the symmetric Cauchy case, if we set a = 1, 

t — we have 



VY^\L]^^\og{2^L) 



1 



TT 



and so we see similar number variance behaviour as in the sine kernel case 1^0.'^) . 

Remark 3. In general, the number variance for the sine kernel determinantal 
process and for Johansson's determinantal process described in the introduction, 
can be expressed in the form H1.8|l with P2t/a°(0j-^) replaced by the functions 

^^^^'^'^i^ and ^^^^'i^i'' + 2a-w'^{cfi+z'^y^ respectively. In the latter case the saturation 
level is again given by the "expectation" term. 

Remark 4. The Cauchy {a = 1) case has the slowest growing non-saturating 
number variance amongst all those considered here. Analogously, the sine ker- 
nel determinantal process has the slowest growing number variance amongst 
all translation invariant determinantal processes whose kernels correspond to 
projections (i.e. the Fourier transform of the kernel is an indicator) see |24| . 

Remark 5. At the other extreme note that as a —> we recover Poisson 
behaviour in that 

liml^r'"m = -(l-e-2^*)- 
a^O a 

1.3 A Poisson approximation for the counting function 

At the beginning of this section we recalled that a system of independent pro- 
cesses (satisfying fairly general conditions) started from a Poisson process on R 
remains in a Poisson process configuration and hence demonstrates a number 
variance linear in L, for all time. Now from 11.6|l and H1.7|l we deduce that 

V"''''''[L]< L/a, y a,c,a,t. (1.13) 

From the integral expression for the number variance 1)1. 3|l , we see that for fixed 
L 

Vj"'^'°[L] ^ L/a ast^oo, 
and for each as L is decreased 



9 



So in both these limitmg cases (as well as in the a ^ case c.f. Remark O 
the maximal linear number variance is attained. The balance of the parameters 
I/, a, a and t, encapsulated by the number variance, determines how "far away" 
from being Poisson the distribution of N"''''"'[0, L] actually is. Below we make 
this observation precise. 

Recall that given a measurable space {il, T) we define the total variation 
distance dTv{'i ■) between two probability measures fii,H2 on ft by 

dTv{f^i,fJ'2) ■■= sup l^i(i^) - ^2(^')|- 

Propostition 1. Let C{N"''^'°'[0, L]) denote the law of the random variable 
7V"''^'°[0, L] defined at H1.4|l . Let Po{L/a) denote the law of a Poisson random 
variable with mean L/a. Then, for L > 1, 

^^^^^^^^<drv{C{Nr^[0, L]), Po{Lla)) < (1 - e'^) ^ " ^^^^ 

Proof. The result is an application of a theorem of Barbour and Hall ^ . Their 
theorem states that if ^ := 1^ is the sum of independent indicator random 
variables indexed by j and qj{L) — E[Ij], A = qj{L) then if we denote the 
law of A by C{A) we have 

-Lmin(l,l/A) < dTv{C[A),Vo{L/ a)) < E'^^W' 

j j 

In our specific case we have N"'''^'°'[0, L] as the sum of independent indicator 
random variables given by H1.4|l . X = L/a and ~ h ~ '^t ''^''°'\^\- D 

Remark 6. For a fixed the Poisson approximation becomes less accurate as 
L —t OQ. The greater the value of a the faster the quality of the approximation 
deteriorates. For a > 1, due to the fact that the number variance saturates, the 
approximation of the law of Af"''^'''[0, L] by a Poisson distribution of mean L/a 
becomes very poor with the total variation distance clearly bounded away from 



2 Gaussian fluctuations of the counting function 

Thus far we have been concerned with the variance of the counting function 
iV"''^'°[0, L] (|1.4() .0f course this variance is, by definition, a description of the 
fluctuation of TV"''^'''^, L] around its mean In this section we will further 
characterize these fluctuations. 

Propostition 2. Let Af"''^'"[0, L], V"''^''^[L] denote, as usual, the counting func- 
tion and number variance. For the cases with a € (0, 1] we have that 

jVr-°[0,L]-L/a 

converges in distribution to a standard normal random variable as L oo. 
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Proof. Recall that the cumulants Ck, fc G N of a real valued random variable Y 
are defined by 



logE[exp(i6iy)] = 



Ck- 



fc=i 



Using the independence of the individual symmetric a-stable processes and then 
applying the Taylor expansion of log(l + x) about zero, we have 

\ogEr[exp{t0N^'''''[O,L])] = ^ log [(e»^ - l)q,{L) + 1 

j=-oo 



m=l 



3 = -oo 



where qj{L) P[X°'''^{t) + Uj e [0, L]] and denotes as usual the underly- 

ing symmetric a-stable process labelled by j with Uj the corresponding starting 
position. Hence, the cumulants of N"'''''^[0, L] are given by 



Ck 



d" ^ 
d^ ^ 



(e*« - 1)^ 



] = -co 



?=0 



(2.2) 



It is straightforward to see that 

CXJ 



give the mean and number variance respectively. More generally from the equa- 
tion H2.2|l it is possible to deduce the following recursive relation 



k—l oo 

771 — 2 j — — 00 



(2.3) 



where I3k,m are constant, finite, combinatorial coefficients which will not be 
needed here. Now let 



iVr'°[0,£]-£/a 



It is easily deduced that the cumulants Cfc, fc e N of Y"''^''^ are given by 



Cl 
Ck 



0, 



Ck 



for fc > 2. 



To prove the Proposition it is sufficient to show that in the limit as i ^ cxd, the 
cumulants correspond those of a Gaussian random variable. That is, we have 
C3 = C4 = C5 = • • • = 0. Equivalently, we need to show 



Ck 



o{{Vt''-'-''[L]f'^) = 0(4/'), for fc > 3. 
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We use an induction argument. 

Suppose that Cm = o{c^^'^) for m = 3, . . . , A; — 1. Assume, without loss of 
generahty, that k is even. We use the inequahty 

fc-1 

q,{L)-q,{L)' = J2q,{Ly-q,{Ly+' 
1=1 

< {k-l)[q,{L)-q,{Lf) (2.4) 

in conjunction with the recursive relation for given by H2.3|l to deduce 
fc-1 fc-1 

^ /3fc,mCm < Cfc < (fc - l)!(fc - 1)C2 + ^ /3fe,mC„. 

m— 2 m— 2 

From our induction supposition this implies that 

0(02^') < Ck < oi4^^) + {k-iy.{k-l)c2. (2.5) 

However, from the results of the previous section, we know that, for these cases 
with a e (0,1], for fc > 3, — V^°'''^'°'[L]^~ — > 00 as L ^ 00. Thus from 

(El Cfe o(c2^^) also. Now using the same arguments as for the inequality 
(EH) we find 

1 < ^ <^ 
^ 3/2 ^ , — ■ 

C2 <=2 JC^ 



Thus we have C3 = o(c2^^). By the induction argument we can deduce that 
(^^"1/2 — > as L — > cx) for all fc > 3 which concludes the proof. □ 

Remark 7. The divergence of the number variance is relied upon in a similar 
way to prove the analogous Gaussian fluctuation results for a large class of 
determinantal processes, see jHlEHll^- In all cases the overall structure of the 
proof is the same. We note that the Proposition could also have been proven 
by applying the Lindberg- Feller Central Limit Theorem (see e.g. |1(J)'). 

Proposition|21 applies to the cases with a E (0, 1]. The following convergence 
in distribution result applies to all cases with a G (0, 2] and is obtained by allow- 
ing both interval length and time t tend to infinity together in an appropriate 
way. 

Propostition 3. For any fixed s E [0, 00) we have that 
iV,"'"'"[0,sti/"]_sti/"/a 



^l/2a 



(2.6) 



converges in distribution as t ^ 00, to a normal random variable with zero mean 
and variance /"'"^'"(s), where 

= r5^^(l-e--("/^)")d.. (2.7) 
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Proof. Since Vf°'''^'°'[st^^"] — > oo as t ^ oo, a similar argument as for the proof 
of Proposition 12 allows us to conclude that 

converges in distribution as t ^ oo to a standard normal random variable. 
Making the change of variable u = st^^"9/a in the integral expression for the 
number variance H1.3|l , using sin^ (u) /v?du = 7r/2 and a double angle formula 
yields 

^^^^^ - ^ r (1 - ^-^^^"/^n - r^^^-w- 



flTT 







Note that from we know /"^^^°(s) < oo for all s < oo and so the result 

follows from the scaling property of the Gaussian distribution. □ 

3 The fluctuation process 

We proceed by adding some dynamics to the fluctuations of the counting func- 
tion and define, for each a £ (0, 2], c > 0, a G R+, the process 

Zt (s) ^175^^ : s e [0,oo). 

3.1 The covariance structure 

We begin to characterize these processes by identifying the covariance structure. 
Lemma 1. {Z"''^'°(s); s e [0,oo)} has covariance structure 

Cov{Zt'^'\r\ ^^'"(s)) = i(r--(s) + r^^\r) - r^^W - s\)). 
Proof. By construction 

N^'^^^iO, {r V s)ii/"] ~ N^^^'^iO, {r A s)t^^"] N^'^^^iO, \r - s\t^/'']. 



ryOL.C.afX W dist r/OL.C.a/ , , \ ryOL.C.a/ . \ 

(k - s|) = [r''^ s)~ Zf - (r A s), 

which implies that 

Var(Z,"'"'°(|r - s|)) =Var(Z4"^"^°(r A s)) + Var(Zf ^"^"(r V s)) 

- 2Cov(Z,"'"'"(r A s), '^'"(r V s)) . 

Rearranging gives 

CoY{zr^'^is),Zr^^{r)) 

= \ (Var(zr'^'''(r A .)) + Var{zr \r V .)) - Var(zr'-'^(|r - .|))) 

On referring back to the definition of /"''^'°(-) we see that this last statement is 
equivalent to the result of the Lemma. □ 

Note that the covariance does not depend on t. 
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3.2 Convergence of finite dimensional distributions 

Given the covariance structure of Z'^''^'"' and the identification of its Gaussian 
one dimensional marginal distributions the natural next step is to consider the 

finite dimensional distributions. 

Propostition 4. Let {G"''^'"(s) : s € [0,oo)} he a centered Gaussian process 
with covariance structure 

Cov{G'''-'-{si),G"'''''{sj)) = ^ir'^^'^Si) + r'^'^'isj) - r'^'%\si - Sj\)). (3.1) 
For any < si < S2 < • • • < < oo we have 

as t OD. 

Proof. As previously noted, the mean and covariance structure of Z^'^'"(s) are 
not dependent on t. Therefore, all that remains is to show that, in the limit 
as t oo, the joint distributions are Gaussian. We again make use of the 
cumulants. 

Recall that given a random vector Y := (Yi, ^2, • • • , Yn) S ffi", the joint cu- 
mulants of Y denoted Cmi.rn2 m„(Y) are defined via the m^'th partial deriva- 
tives of the logarithm of the eharacteristie function of Y. That is, 

f d \ "/ d \ / 9 \ 



If 



^OA...,^,...,o(Y) = E[y,]=o 

i'th 

Co,o,...,^,...,o(Y) = Var[y,] = E,, 

i'th 

<^o,....^,...,o,...,^,...,o(Y) = Cov[y„y,] = s.,- 

i'th j'th 

and in particular 

n 

(Y) = 0, whenever ^ mj > 3, 

then Y has a multivariate normal(0, S) distribution. To prove the Proposition 
it is enough to show that 

'"'"(sa)) ^ MultivariateNormal(0, S"''^'") (3.2) 

in distribution as t ^ oo, where 1]"''=''^ is the 2x2 covariance matrix 

( /"'^'"(si) M/"''''tsi)+r'"'"(s2)-r'^'isi-s2^)A 

U(/"'''°('Si)+r'''"N-r'^'isi-s2^) /"■^■"(s2) )' 
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We begin by computing the characteristic function of 

Using the independence of the individual particles we have 
Ep [exp(i0iiVf '"'"[0, si^i/"] + i^aiVf '"'"[0, ssi'/"]) 
= Ep [exp + e2)Nt''''[0, siii/"] + 02Afr'''"[sii'/", S2i'/" 

oo 
j=-oo 

For ease of notation we will henceforth let 

qjisi,Sr) F[Xf'^{t)+U, e < S/ < < C50. 

The joint cumulants are given by 



9 \ ™2 



log 



j=oo 

Now using the fact that 

d \ '"v d \ ™2 



e'(^i+«2)gj(0,si) VTOi,m2 s.t. mi > 1 



and 



^) e^'^+'^^kjiO, ^1) + e''^g,(.i, .2) + 1 - 9,(0, .2) 



along with 



= 1, 



0=0 



we deduce the following generalization of the recursive relation H2.3|l . with ob- 
vious short-hand notation for the joint cumulants 



C*7m,m2 — ^ ^ Pk.l,nii .m2 ^k,l 

kd: 

2<fe+i<mi+m2-l 

00 

+ (mi +7712-1)! 9/0,si)-(7,(0,si)"%(0,S2)™^ 
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Now suppose that Ck,i — o{t ) for all fc, I such that k + I e {3, 4, ... , toi + 
m2 — 1} and without loss of generality assume that mi + m2 is even. Since 

< 9,(0,.si) -9,(0,si)'"^rz, (0,S2n < 9,(0, si), 
if we assume that the above induction hypothesis holds, then we have 

oit""^^^) < C„„,„, < o{t''^^^) + (mi + m2 - (3.3) 

a 

which implies that Cmi,m2 = o(t as ) also. We check the third order joint 
cumulants directly and deduce that 

o{t^) < Ck,i< o{t*^) + 2[-^——j whenfc + ; = 3, 

since the variances and covariance of {N^-°-[0,sit^/°'],N^-°-[0,S2t^^'^]) (i.e. 
(^2,0, Co,2, Ci^i) grow at most like t^/" as t ^ oo. Therefore, by induction, 
whenever mi + m2 > 3 we have 

7 ; 7-T- > as J ^ CX). 

In terms of the joint cumulants of z"'''''^ this implies 

Cmi,m2iZ^'^'°'{si), Z^''^'"'{s2)) ^ as i ^ oo whenever mi + TO2 > 3, 
from which the claim H3.2|l and the Proposition follow. □ 

3.3 A functional limit for the fluctuation process 

In order to give a functional limit result we consider a continuous approximation 
{Z"''^''^{s) : s e [0,cxd)} to the process {Z"''^'"(s) : s e [0, oo)}. Let 

Zr'^i^) ■■= ^ ^"''^1/2^ (3-4) 

where N"''^''^[0, st^^"] is defined to be equal to N^''^'"'[0, st^^°'] except at the 
points of discontinuity where we linearly interpolate. Let C[0, 1] be the space of 
continuous real valued functions on [0, 1] equipped with the uniform topology. 
We shall denote the measure induced by {Z"'^'°(s) : s € [0, 1]} on the space 
(C[0, 1], S(C[0, 1])) by Q"''^'". To simplify notation we restrict attention to 
the interval [0, 1] but note that the ensuing functional limit theorem extends 
trivially to any finite real indexing set. The remainder of this section is devoted 
to establishing the following weak convergence result. 

Theorem 4. Let Q°' '^-°- be the law of the centered Gaussian process {G"''^'"(s) : 
s G [0,1]} introduced in the statement of Proposition^ Then 



Qa,c,a 
t 



as t ^ OO. 
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Proof. Note that by definition 

izr'"(s)~zr''"(s)i<^. (3.5) 

Thus, as t ^ oo, the finite dimensional distributions of Z"''^'°(s) must converge 
to the finite dimensional distributions of the limiting process G'"''^'°(s) to which 
those of ^"'"^'^ converge. Hence, immediately from Proposition 0] we have that 
for any < Si < S2 < ' ' ' < Sn < oo, 

as t — > cx). Therefore, by a well known result of Prohorov's (see e.g. jHI), the 
proposed functional limit theorem holds if the sequence of measures {Q"'*^'"} 
is tight. Indeed this tightness requirement follows from Proposition [S] given 
below. □ 

The sufficient conditions for tightness verified below are stated in terms of 
the distributions 

P(Zf eA)^ g^'"'"(A) A e 6(C[0, 1]), 

and the modulus of continuity, which in this case is defined by 

^zr■^<5):- sup izr''^(.)-zr'"(r)i, <5e(o,i]. 

\s-r\<S 

Propostition 5. Given e, A > 3(5 > 0, to E N such that 
P[w(Z4"'"'°,(5) > A] < e, fort>to. 
Proposition 13 is proven via the following scries of Lemmas. 
Lemma 2. Suppose 0<u<r<s<v<l, then 

|Zf '"^°(s) - '"^°(r)| < |Zf '"'"(w) - '"'"(u)! + {v- u)t^^^". 
Proof. Clearly, by construction 

< 7Vt"'"'"[0,sii/"] - iVf'"'"[0,rii/"] < 7Vf'"'"[0,wti/"] - iVf '"'"[0, ut^/"]. 
Therefore, using the definition of Z"''^'", we have 

a a 

The result follows by rearranging, using the facts a G R+, (w — u) > (s — r) > 
and then considering separately each case 

Zj"^"'"(s) - Zj"^"^°(r) > or Zj"^"^°(s) - '"^"(r) < 0. 

□ 

Lemma 3. 

\Zr^is) - Zr^^ir)\ + l^r'""(=s) - Zr''^ir)\. (3.6) 
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Proof. Follows from H3.5|l and an application of the triangle inequality. □ 

To obtain results on the distribution of the modulus of continuity for our 
sequence of processes we divide the interval [0, 1] into m disjoint subin- 

tervals of length approximately S as follows. Let 

= ro < rfc, < • • • < rfc„_i < rfc„ = 1, (3.7) 

where we define 

k, := j\Sf/''], je{0,l,2,...,m-l} 
and [•] denotes the ceiling function. We have 

5<rk^-rk,_, <S+^ J e {l,2,...,m-l}, 

with the subintervals [ri_i,ri], i e N typically being shorter in length. 
Lemma 4. 

m— 1 , „ 

Proof. Given the partition H3.7|l . standard methods (see Theorem 7.4 of [HI) 
yield 

m— 1 

p\w{Zr'\S)>x\<J2p\ sup \Zr-''{s)-Zr-''{r,^)\>^]. (S.S) 

■ n rk ■ <s<rk . , , "J -' 

By the triangle inequality we have 

\zr^%s) - zr^%r^^)\ < izr-'is) - zr^''in)\ + izr^n) - zr^''ir,^)\. 

(3.9) 

Now if s G [j'/cj , ^'fej+i], then either 

''^'"(s) = ''='"(r,) for some i G N 
immediately simplifying H3.9|l . or 3i G N such that 

Tk, < n-1 < s < r^ < rkj^i 
in which case from Lemma 12 we have 

izr'°(s)-zr'''(n)i 

< |Zr^^"(rO - Zr'^''^(r,_i)| + (r, - r,_i)ti/2« 

= \zr-''in)-zr^^{r,^,)\ + ^ 
\zr'\n) - zr-\rk,)\ + i^r"^°(^.-i) - ^r'^''^(^fe.)i + (3.10) 



< 
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Therefore, using the inequahty (|3.9|l m conjunction with H3.1Uf) and Lemma 13 
we have that for s, e [rk^ , rkj^^] 



< 



< 3 max 



< 3 



ryOL.C.a ( \ ryOL.C.a / \ 



+ 



3tl/2o 



Thus 



sup 



< 



max 

J < t < J -I- 1 



> 



7 A 

> 



3ti/2a - 9 J 

Substituting this last inequahty into H3.8|l gives the statement of the Lemma. □ 

Now that we have reduced the study of the distribution of the modulus of 
continuity to that of the maximum fluctuation over our constructed subintervals 
we can progress by introducing a maximal inequality. In order to do this we use 
the following known result taken from 6 and paraphrased for use here. 

Theorem 5. Consider a sequence of random variables {S,q}q>i and the associ- 
ated partial sums 

u 

Su ■— ^ S,q So :— 0. 

g=l 



Let 



If 



M, 



n .- max 15*^ 

l<«<r! 



Q < u < V < n 



u<l<v 



for 7 > and some 6i, 62, ■ ■ • , G then 

' 0<l<n 

where k is a constant. 
Proof See 0, Theorem 10.2 
Lemma 5. 

P[ maxjZr-''(r,^.+„) - Zr-^{r,^)\ > 7] < ^{nt~'^"'y 
where k is constant and 7 > 0. 



□ 
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Proof. Let 
Then 

U 

Su ■=^£.q = Zt-'^-^irk^+u) - Z^'^'^'irk^) So 0. 

q=l 

In this case, for v,u G N, applying Chebyshev's inequaUty and using the defini- 
tion of /"''^'°(-) and the upper bound given at H1.13|l . we have 

P[\S. - Su\ > 7] = P[|Zr^'"(rfe,+„) - Zr'''{rk,+u)\ > 7] 

7^ 

^ (t)-u)i-i/" 
~ a7^ 

— 9 

a7^ 

Thus we can take bi — - — ^ — for I = I, 2, . . . ,n and apply Theorem which 
gives the maximal inequality of the Lemma. □ 

Concluding the proof of Proposition is now straightforward. 
Proof of Proposition\^ Taking n — fcj+i — kj in the statement of Lemma|Slgives 

07^ 

Substituting this last inequality with 7 = (■§■ ^ 3^72^) (which is strictly positive 
for sufficiently large t) into the inequality given by Lemma 01 gives 

,[„,(zr".^)>A]<g^'^''^'°"'"''°/. 

On solving the appropriate quadratic equation we find that we can make this 
upper bound less than e by choosing 6 from 



_ ^-1/2q 



(A - 21t-l/2a) ^ -t-l/2a + (A - 21i-l/2")) fl (0, 1) , 



where Cg = a/ gimK • Since this interval is non-empty for sufficiently large t this 
completes the proof. □ 

3.4 Properties of the limiting process G°''^'"'{s) 

We have constructed a family {{G°'''''°'{s), s e [0,cx))),q; e (0, 2],c > 0, a G 
M+} of centered, continuous, real- valued Gaussian processes with the inherited 
covariance structure 

CoviC^'^^'^is), C^-^'^ir)) = i(r^^^''(s) + /"•^■''(r) - /"•^■''(Is - r|)) . (3.11) 
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It is clear that the processes are recurrent for 1 < a < 2 as the number variance 
saturates and the stationary distribution is normal with mean zero and variance 
Ksat{c(, c, a, 1). We are able to deduce further properties of these limit processes 
by using our earlier results on the number variance for the systems of symmetric 
a-stable processes from which they are constructed. 

Propostition 6. G"''^'"(s) is non-Markovian. 

Proof. Recall (see e.g. ^H], Chapter 13) that a Gaussian process with indexing 
set T C R and covariance function p : T i-^ M. is Markov if and only if 

P{s,u)^ ^ ys,u,r(^T. 

p{r,r) 

It is clear that this relationship does not hold in general for the covariance 
function lprTT|l . □ 

From the results of previous sections, since /"^'^'"(s) = V^]"''^'''[s], we know 
that 

ja,c,a/ /^_g-2cx£ as a ^ 

a 

and 

/ ' ' (s) ^ ka c a S. 

Therefore, G"''^^" appears to start out, for small "time" as a scaled Brownian 
motion and as a ^ this initial Brownian character prevails for longer. We 
capture this more precisely in the following easily verified proposition. 

Propostition 7. 1. {G"''^'°(s) : s e [0,1]} converges weakly to a Brownian 
motion {B(i^^s) : s e [0, 1]} as a ^ 0. 

2. Let G^''='"(s) = e-i/2G"^'='°(es). Then {Gf'=^°(s) : s G [0,1]} converges 
weakly to a Brownian motion {B{^) : s G [0, 1]} as e 0. 

Remark 8. The covariance structure of G"''^'" is similar to that of a Brownian 
bridge. Recall that the standard Brownian bridge {B''^{s),s € [0, a]) of length 
a, is a centered Gaussian process with covariance structure 

Cov{B'"'{s),B''''{r)) = sAr - ^ 

and arises as a weak limit of many empirical processes. In particular, it may be 
obtained from the appropriately scaled counting functions of a Poisson process 
on M (see e.g. ^Hl)- We can re- write the covariance (|3.1|l in the alternative form 

CoviG'^-^^-^{s),G"'^''^{r)) = - /° /%2/a»(y,^)dydz, 

a Jo Jo 

but we see that a precise match would require P2/a° (j/i z) = 1. 

Propostition 8. The process G"''^'" has stationary increments which are neg- 
atively correlated. 
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Proof. It is straightforward to see that the increments have zero mean and that 
for any s, r G [0, oo) 

Var(G'"''='"(s) - G"'"'°(r)) = /"'"'""{Is - r|). 

In addition, for u > we have 

Cov(G"^"^°(s) - C'^'^iO), G"''=^°(r + s + u) ~ G"'"'"(s + u)) 

= I {/"•'•"{s + r + u)- r-^-^ir + u)- [r^=^"(,s + u) ~ /"^^^''H]) . (3.12) 

Since /">='° is a concave function 

/"■"■'"{s + r + u)- /"'^■°(r + u) < /"''='"(s + u) - /"■'^'''(u), 

for all s,r, u e [0, oo) so it follows that this covariance is non-positive. □ 

Propostition 9. G"''^'" is not in general self- similar. For any constant 6 G M 
we have the relationship 

C'^'^ibs) b^^^G^-'T^-^is). 

Proof. Both sides of the proposed equation have zero mean and a Gaussian 
distribution. It is clear from the expression given for f"''^''^ at 1)2. 7|l that the 
variances/covariances also agree. □ 

Propostition 10. G"''^'° is a long memory (or long range dependent) process 
in the sense that the covariance between increments decays as a power law as 
the separation between them is increased. More precisely, for a £ (0, 2) we have 

Cov{G"^''^''is) ~ G"^"^°(0), G"''='"(r + s + u) - G"^'=^°(s + u)) ku-^°'+^\ 
where k is a constant depending on a, c, a, s and r. 

Proof. The covariance in question is expressed in terms of the function 
at H3.12(l . Note from l|1.9|l that we already know the asymptotic behaviour of 
the individual components of this expression. Applying I'Hopital's rule twice in 
succession yields the given power law. □ 

We have already mentioned a similarity between the covariance structure 
of G"''^'° and that of Brownian motion. More generally we can draw parallels 
between our limiting process and fractional Brownian motion. Recall (see for 
example ^Hl) that a fractional Brownian motion {Wh{s),s > 0) with Hurst 
parameter H £ (0, 1) is a centered, self-similar Gaussian process with covariance 
function 

Cov{Wh{s), Wnir)) - ^ [s'" + r^" - \s - r\^") , (3.13) 

The case H — ^ corresponds to a standard Brownian motion. 

Note the resemblance between the form of the covariance functions (|3.11|l 
and H3.13|l . Heuristically, we can deduce that /"■'^•"(s) may be approximated 
by a function of the form Ka,c,aS^"^°'''''" where Ha^c,a '■ [0, oo) i— > [0, |] is 
a monotonically decreasing function with initial value Ha^c,a{0) — ^- Thus 
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loosely speaking G°' '^ °' can be viewed as a type of fractional Brownian motion 
with time varying Hurst parameter. In particular, the long range dependence 
property of Proposition 1101 may be compared to the analogous statement for 
fractional Brownian motion: 



Cov{Wh{s) - Tyif(O), Wnir + s + u)-Wh{s + u)) ku^"-"^, 

where A: is a constant depending on H, s and r. We make the link between 
the process G"'*^'" and fractional Brownian motion precise with the following 
statement. 

Propostition 11. For a S (0, 1) let 

Then 

{G"''='''(s) : s e [0, c5o)} ^ {kUl^ Wi^ (s) : s e [0, oo)} as h ^ oo 
where ka^c,a = |f r(a - 1) sin ( - a7r/2) . 

Proof. From (|2.7|l . by applying a Taylor expansion we deduce that 

61-" ' air Jq 

.... „i«-2 

Now by the Dominated Convergence Theorem, 

+ 1 .... 7,J"-2 

3=2 



■/ sin^K2)V ^ \ (2c)^V-^^"^ du. 

J — ^ 

Convergence Theorem, 
h^oo air Jq ^ ] ! oJ" " 



hb(u) 

^ ± / sin2(u/2)5] lim L^(2c)^V-2,a« 
= 

since, setting M{u) — siv?{u/2){l — exp(— 2c(m/s)"))/u^ + siv?{u/2) /v?'"" , we 
have a positive integrable function such that |/i6(u)| < M{u) for all 6 G R. This 
implies that 



lim — = ka.c.a s 



which allows us to conclude 



6- 



k 



s — r 



2 
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The processes are Gaussian therefore the convergence of finite dimensional dis- 
tributions is implied by the convergence of the covariance functions and tightness 
follows easily from, for example, |16| Corollary 16.9 and well known expressions 
for the even moments. □ 

Remark 9. We mention that various other long-range dependent Gaussian 
processes have recently been found to arise from the fluctuations of spatially 
distributed particle systems, see JT] and references within. Notably, in this 
context, the spatial particle configurations of infinite systems of symmetric a- 
stable processes started from a Poisson process on M have been considered. In 
these cases, fractional Brownian motion with Hurst parameter iJ = 2 — i, a € 
(1,2] was obtained as a scaling limit of the occupation time process (essentially 
scaling the counting function in time rather than in time and space as in this 
paper). 

Remark 10. It seems natural to ask whether, in the same fashion as we created 
Qa,c,a ^ similar limiting processes could be constructed from Johansson's systems 
of non-colliding Brownian motions. Unfortunately, the formula for the averaged 
number variance given for these processes in 1141 does not scale in time and 
space in the same convenient way as 1^"''^'''[L] in this case. However, as noted by 
Johansson, letting i — > oo in his model, one obtains the sine kernel determinantal 
process, from which a limiting Gaussian process (parameterized and scaled in a 
completely different way) was constructed in |23| . 
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